Non-local interactions in hydrodynamic turbulence at high Reynolds numbers: 

the slow emergence of scaling laws 
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We analyze the data stemming from a forced incompressible hydrodynamic simulation on a grid 
of 2048^ regularly spaced points, with a Taylor Reynolds number of R}, ~ 1300. The forcing is given 
by the Taylor-Green flow, which shares similarities with the flow in several laboratory experiments, 
and the computation is run for ten turnover times in the turbulent steady state. At this Reynolds 
number the anisotropic large scale flow pattern, the inertial range, the bottleneck, and the dissipative 
range are clearly visible, thus providing a good test case for the study of turbulence as it appears 
in nature. Triadic interactions, the locality of energy fluxes, and structure functions of the velocity 
increments are computed. A comparison with runs at lower Reynolds numbers is performed, and 
shows the emergence of scaling laws for the relative amplitude of local and non-local interactions in 
spectral space. The scalings of the Kolmogorov constant, and of skewness and flatness of velocity 
increments, performed as well and are consistent with previous experimental results. Furthermore, 
the accumulation of energy in the small-scales associated with the bottleneck seems to occur on a 
span of wavenumbers that is independent of the Reynolds number, possibly ruling out an inertial 
range explanation for it. Finally, intermittency exponents seem to depart from standard models at 
high Rx, leaving the interpretation of intermittency an open problem. 

PACS numbers: 47.27.ck; 47.27.Ak; 47.27.Jv; 47.27.Gs 



I. INTRODUCTION 



Turbulence prevails in the universe, and its multi-scale 
properties affect the global dynamics of geophysical and 
astrophysical flows at large scale, e.g. through a non-zero 
energy dissipation even at very high Reynolds number 
Re. Furthermore, small-scale strong intermittent events, 
such as the emergence of tornadoes and hurricanes in at- 
mospheric flows, may be very disruptive to the global 
dynamics and to the structure of turbulent flows. Typ- 
ically energy is supplied to the flows in the large scales, 
e.g., by a large scale instability. The flow at these scales is 
inhomogeneous and anisotropic. In the standard picture 
of turbulence, the energy cascades to smaller scales due 
to the stretching of vortices by interactions with similar 
size eddies. It is then believed that at sufflciently small 
scales the statistics of the flow are independent of the ex- 
act forcing mechanism, and as a result, its properties are 
universal. For this reason, typical investigations of turbu- 
lence consider flows that arc forced in the large scales by 
a random statistically isotropic and homogeneous body 
force P, 01- However, how fast (and for which measured 
quantities) is isotropy, homogeneity, and universality ob- 
tained is still an open question. 

The return to isotropy has been investigated thor- 
oughly in the past, by analysis of data from experiments 
and direct numerical simulations (DNS) [1, 0, H, d, 0; 111 ■ 
However, lack of computational power limited the numer- 
ical investigations of anisotropic forced flows to moderate 
Reynolds numbers, for which a clear distinction of the in- 
ertial range from the bottleneck, and from the dissipative 



range, cannot be made. Only recently the fast increase 
of computational power permitted DNS to resolve suffl- 
ciently small scales, such that a flow due to an inhomo- 
geneous and anisotropic forcing develops a clear inertial 
range with constant energy flux. As a result, this kind 
of questions can be addressed anew. To give an estimate 
of the size of the desired grid, we mention that in recent 
simulations Q an incipient inertial range was achieved 
for a resolution of 1024'^ grid points, while for a 512'^ run 
the range of scales between the large scale forcing and the 
bottleneck was much less than an order of magnitude. In 
all cases, the flow was resolved since kmax/kri ^ 1, with 
kmax the maximum wavenumber in the simulation and 
k^i the dissipation wavenumber built on the Kolmogorov 
phenomenology. 

Of particular interest in the study of turbulent flows 
is the issue of universality It is now known that two- 
dimensional turbulence possesses classes of universality 
[l0| . and at least for linear systems such as the advec- 
tion of a passive tracer, there is evidence of universality 
of the scaling exponents of the fluctuations How- 
ever, recent numerical simulations of three dimensional 
turbulence Q showed that scaling exponents of two dif- 
ferent flows (one non-helical, the other fully helical) were 
measurably different at similar Reynolds number. It is 
yet unclear whether this is an effect of anisotropics in 
the flow, or of a finite Reynolds numbers. If this is a 
finite Reynolds number effect, one then needs to ask how 
fast its convergence to the universal value is obtained. 
If the convergence rate is sufficiently slow then finite 
Reynolds effects should be considered when studying tur- 
bulent flows that appear in nature, at very large but flnite 
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Reynolds numbers. Thus, the question of the universal 
properties of turbulent flows at high Reynolds numbers 
remains somewhat open. 

The recovery of isotropy, the differences observed in 
the scaling exponents, and the slow emergence of scaling 
laws have been recently considered in the context of the 
influence of the large scales on the properties of turbu- 
lent fluctuations dZ, dM ■ The study of nonlocal inter- 
actions between large and small scales has been carried 
in experiments and in simulations 0, [l^, [ij, [l^, [3, [13, 
in, [IS H EH M, M, IM, M, M, M at small and mod- 
erate Reynolds numbers. In simulations with 1024'^ grid 
points [i^l, it was found that although most of the flux 
is due to local interactions, non-local interactions with 
the large scale flow are responsible for « 20% of the total 
flux. It is however unclear how the amplitude of these 
interactions scale with the Reynolds number. 

In this context, we solve numerically the equations for 
an incompressible fluid with constant mass density. The 
Navier-Stokes equation reads 



atu 4- u • Vu = -VP 



(1) 



with V • u = 0, where u is the velocity field, V is the 
pressure divided by the mass density, and v is the kine- 
matic viscosity. Here, F is an external force that drives 
the turbulence. The mode with the largest wavevcctor 
in the Fourier transform of F is defined as kp, with the 
forcing scale given by 2TT/kp. We also define the viscous 
dissipation wavcnumbcr as kri = {t/v^Y^^ , where e is the 
energy injection rate (as a result, the Kolmogorov scale 
is 77 = 27r/fc^). 

The results in the following sections stem from the 
analysis of a series of DNS of Eq. ^ using a parallel 
pseudospectral code in a three dimensional box of size 
27r with periodic boundary conditions, up to a resolution 
of 2048^ grid points. The equations are evolved in time 
using a second order Runge-Kutta method, and the code 
uses the 2/3-rule for dealiasing. As a result, the maxi- 
mum wavenumber is kmax = N /Z where N is the number 
of grid points in each direction. 

With L and A defined as 



27r 



jE{k)k-^dk 
j E(k)dk ' 



A = 27r 



jE{k)dk 
J E{k)k^dk 



1/2 



(2) 



the integral scale and Taylor scale respectively, the 

Reynolds number is Re = UL/v and the Taylor based 

1/2 

Reynolds number is R\ = U\/v. Here, U = {v?) is 
the r.m.s. velocity and E{k) the energy spectrum. The 
large scale turnover time is T U /L. Note that, with 
these definitions. Re and Rx used in this paper are larger 
than the ones stemming from the definitions used by the 
experimental community (see e.g., [1^) by a factor of 
2^(3/5)1/2 7» 4.87. 

Simulations were done with the same external forcing 
(see Tabic U for the parameters of all the runs), with 
U ~1 \n all steady states. The forcing F corresponds to 



TABLE I: Parameters used in the simulations. A'^ is the linear 
grid resolution, v the kinematic viscosity. Re the Reynolds 
number, and R\ the Taylor based Reynolds number. 



Run 


iV 


V 


R. 


Rx 


I 


256 


2 X 10"^ 


675 


300 


II 


512 


1.5 X 10"^ 


875 


350 


III 


1024 


3 X 10"'' 


3950 


800 


IV 


2048 


1.2 X 10"* 


9970 


1300 
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FIG. 1: (a) Energy spectrum in run IV compensated by k '^^^ . 
The inset shows the energy flux, (b) Energy spectrum in runs 
III (dotted) and IV (sohd) compensated by fc"*^"^. Wavenum- 
bers are normalized by the dissipation wavenumber fc,,. 



, Taylor-Green (TG) flow 

F = /o [sm{kpx) cosikpy) cos{kpz)x— 
— cos{kpx) sin(kpy) cos{kpz)y] , 



(3) 



where /o is the forcing amplitude, and kp = 2. The 
turbulent flow that results has no net helicity, although 
local regions with strong positive and negative helicity 
develop. 



II. THE SLOW EMERGENCE OF A 
KOLMOGOROV-LIKE SCALING 

We first concentrate on the global dynamics of the 
2048"^ run (run IV). Figure [IJa) shows the compensated 
energy spectrum in this run, as well as the corresponding 
energy fiux Il{k), both taken in the turbulent steady state 
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after the initial transient. The energy flux is constant 
in a wide range of scales, as expected in a Kolmogorov 
cascade, but the compensated spectrum has a more com- 
plex structure in that same range of scales. The salient 
features of this spectrum are well-known from previous 
studies. Small scales before the dissipative range show 
the so-called bottleneck effect with a slope shallower than 
^-5/3 _ ^Yic other hand, larger scales have a spectrum 
with a slope slightly steeper than k~^^^, an effect that 
is even clearer in the simulation performed at larger spa- 
tial resolution [l| on a grid of 4096'^ points; this small 
discrepancy with a Kolmogorov spectrum is attributed 
to intermittency, i.e. to the spatial scarcity of strong 
events leading to non-Gaussian wings in the probability 
distribution functions of velocity gradients. 

The bottleneck effect is not fully understood but 
clearly corresponds to an accumulation of energy at the 
onset of the dissipation range. It has been attributed to 
the quenc hing of local interactions close to the dissipa- 
tive scales dfl Hll, [H, , or to a cascade of helicity [33| 
whose energy spectrum would follow a k~^/^ power law. 
The quenching of local interactions in the bottleneck was 
measured directly in simulations in and will be also 
shown here for run IV (see below. Figs. [2]|4]). The k~"^^^ 
spectrum is also compatible with the present data, as 
shown in Fig. [ijb) giving the energy spectra in runs III 
and IV compensated by k~~^^^. However, we observe that 
the width of the bottleneck appears to be independent of 
the Reynolds number; this indicates that the origin of 
the bottleneck is more likely a dissipative viscous effect 
than an inertial range effect. If helicity plays a role in 
the formation of the bottleneck, it has to be connected to 
the local generation of helicity at small scales due to the 
viscous term in the Navier-Stokes equation. Purely he- 
lical structures are exact solutions of the Navier-Stokes 
equation, and as a result an increase of helicity in the 
small scales could quench local interactions and the cas- 
cade rate (as assumed in Ref. [3^]). 

The relative strength of local versus non-local interac- 
tions between Fourier modes in the shell-to-shell transfer, 
and in the energy flux can be measured in numerical sim- 
ulations with the help of a variety of transfer functions 
[H, m, [li, [13, ill. Specifically, the amplitude of the 
basic triadic interactions between the modes in shells 
P and Q is defined as: 

T^{K, P,Q)^-JuK- (up • V)uQdx3, (4) 

where the notation uk denotes the velocity field fil- 
tered to preserve only the modes in Fourier space with 
wavenumbers in the interval [K,K + 1). Picking a 
wavenumber in the inertia! range (here Q = 40), we show 
in Fig. [2Ia) its amplitude as a function of P and K — Q 
for run IV. Specific values of two levels are indicated 
as a reference (the maximum, indicated by the arrow, 
corresponds to P = kp)- As a comparison, in run III, 
max{T3(ii:,P,Q = 40)} sa 1.4 x lO'^ indicating that a 
decrease of the relative amplitude of the non-local triadic 
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FIG. 2: (a) Amplitude of the triadic interactions Ts{K , P, Q) 
for Q = 40 as a function of K and P in run IV. (b) Sliell-to- 
shell energy transfer Tz{K, Q) in the same run; several values 
of Q are superimposed. 



interactions with the large scale flow {P = kp) occurs as 
the Reynolds number increases. However, the non-local 
coupling of the modes with P fa kp is still dominant in 
run IV. 

The relevance of these interactions in the transfer of 
energy between scales can be quantifled by studying the 
shell-to-shell transfer and the net and partial fluxes. The 
energy transfer from the shell Q to the shell K, integrat- 
ing over the intermediate wavenumber, is defined as: 

T^iK, Q) = Y. n{K, P,Q)^- f UK-iu- V)uQdx3 . 

p 

(5) 

It has the same qualitative behavior as in runs at lower 
Reynolds number [see Fig. [HJb)]. The minimum of T2 for 
K ^ Q M —kp for all values of Q, and the maximum for 
K — Q « kp, both denote that the energy is transfcrcd 
from the nearby shell K—kf to the Q shell, and transfered 
from this shell to the nearby shell K + kf. As a result, 
as we increase the Reynolds number, the shell-to-shell 
energy transfer is still local but not self-similar, mediated 
by strong non-local triadic interactions with the large 
scale flow at fcp d, [H, M, M, M, 113, HI, HI . 

It has been observed that although the individual 
non-local triadic interactions are strong, as modes are 
summed to obtain the energy flux, non-local effects be- 
come less relevant. To quantify further the net contribu- 
tion of the local and non-local effects to the energy flux. 
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FIG. 3: (a) Ratio of large-scale to total energy flux 
Tl-Ls{k) /Tl{k) as a function of wavenumber for runs I (dash- 
dot), II (dash), III (dot), and IV (solid), (b) Ratio of non- 
local to total energy flux nNL(fc)/II(A;) for the same runs; 
wavenumbers are in units of fc,,. The fluxes are deflned in 
Eqs. ((6}, 0, and JH). 

we introduce the total flux 

k 

nW = -EEE^3(if,/',Q), (6) 

K=0 P Q 

the energy flux due to the non-local interactions with 
only the large scale flow 

k 6 

nL5(fc) = - E E E P, Q), (7) 

K=OP=Q Q 

and the energy flux due to all the interactions outside the 
octave around wavenumber k (i.e., all non-local interac- 
tions) 

k k/2 

nNL(fc) = -^ 5^^T3(if,P,Q). (8) 

A'=OP=0 Q 

Figure[2ja) shows the ratio nL5(fc)/n(fc) as a function 
of wavenumber for run IV. The same ratio for the lower 
resolution runs in Table U are also shown here as a ref- 
erence. If the cascade is due to local interactions, this 
ratio should decrease as smaller scales are reached. We 
observe however that, at small scales, a plateau obtains 
within which this ratio remains relatively constant. This 
is observed in runs III and IV, the two runs at the highest 



Reynolds numbers. Note also that the plateau lengthens 
as Rx increases: the length of the plateau corresponds 
roughly to the length of the inertial range (including the 
bottleneck) at those Reynolds numbers. Finally, the am- 
plitude of the plateau decreases as the Reynolds number 
is increased, indicating a smaller contribution of the in- 
teractions with the large scale flow, relative to the total 
flux. A detailed study of its dependence with Reynolds 
number is discussed in the next section. 

As previously mentioned, the ratio Hls/II does not in- 
crease in the range of wavenumbers that spans the bot- 
tleneck. It is the contribution of all non-local interac- 
tions (interactions with all the modes outside the octave 
around a given wavenumber k) that becomes dominant 
in this range. Figure [SJb) shows the ratio Il-^L{k)/Il{k) 
for the runs in Table[T](note the wavenumbers are plotted 
in units of fc^). As scales closer to the dissipativc range 
are considered, the contribution of all the non-local inter- 
actions increases, in agreement with the findings in Ref. 
jsot - Moreover, the amphtude of Uj<!L{k)/Il{k) when k 
is in units of fc^ is roughly independent of the Reynolds, 
in agreement with a width of the bottleneck independent 
of Re and controlled by the growth of nNL(fc) as k gets 
closer to the dissipation scale. 

It is worth noting that even at the highest Reynolds 
number examined here, there is still a significant con- 
tribution of nonlocal interactions (IIls and IlNi) to the 
total energy flux in the inertial range. The comparison 
with runs at smaller resolution shows a qualitative agree- 
ment and the persistence of the described non-local ef- 
fects. What the new computation at Rx ^ 1300 allows, 
though, is to determine the scaling of the relative im- 
portance of nonlocal effects in Navier-Stokcs turbulence 
when the Reynolds number is increased, as we discuss 
next. 



III. SCALING LAWS IN TURBULENT FLOWS 

Numerical simulations do not excel in the determina- 
tion of scaling laws in turbulent flows. The resolutions 
allowed by present day computers barely allow for the 
existence of a well-defined inertial range. Indeed, the ob- 
servation of Fig. [1] shows that, at this Taylor Reynolds 
number, the Kolmogorov inertial range covers less than 
one order of magnitude in scale (although, as noted be- 
fore, the fiux is constant in a larger range of scales). This 
could be an indicator that solutions more complex than 
simple power laws hold in the inertial range (40| . The pio- 
neering computations of the Japanese group on the Earth 
Simulator using random forcing has allowed, however, for 
some scaling laws to emerge, although, as these authors 
observed, not all physical quantities of interest converge 
to asymptotic values at the same rate [4l|, . We here 
display such scaling laws for the particular flow studied, 
namely the Taylor-Green flow, relevant to several labora- 
tory experiments. In particular, we are interested in the 
scaling of the relative amplitude of local and non-local 
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FIG. 4: Scaling of (a) the flux ratio IlLs/n and (b) the non- 
local flux ratio IlNL/n as a function of Reynolds number. 
Both ratios are evaluated at the Taylor scale, and several 
slopes are indicated as references. 
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FIG. 5: Kolmogorov constant C as a function of Reynolds 
number; the solid line gives the best fit C = 4.60i?;:'' " +0.64. 



interactions, as well as other quantities often studied in 
the context of turbulent flows, whose sca ling will be used 
as a criteria to classify the runs [U, [13, l43| |. It is worth 
mentioning in this context that, with four runs, we can 
only show the results are consistent (or at least, not in- 
consistent) with a particular scaling. 

Figure Hl^a) gives the scaling of the flux ratio 
nLs(fc)/n(fc) with the Reynolds number. To this end, 
we take the Taylor scale A as a reference scale in the in- 
ertial range, and we evaluate Ili^s{k)/Il{k) for each run 
at the Taylor wavenumber kx = 27r/A. The best fit to all 
the runs gives XIls/II ~ -Rj^ '', although the dependence 
of the ratio Hls/II with Re seems to change slightly for 



FIG. 6: (a) Kurtosis and (b) skewness of the velocity incre- 
ments as a function of the Taylor based Reynolds number. 
Results are given for two different increments: the Taylor 
scale and the dissipation scale (*). The slopes indicated 
as a reference are from experimental results. 



Run IV. A best fit of the last two points (runs III and IV) 
gives a dependence ~ R^^'^ (as was discussed in Sect. |IT1 
these two runs show a developed inertial range). 

We also evaluate the ratio IlNi/n at the Taylor 
wavenumber; its dependence with the Reynolds number 
is shown in Fig. HJ^b). Here, the ratio in runs III and IV 
is compatible with a slower decay ~ R~^ '^. The anoma- 
lous behavior of runs I and II in Fig. [DJb) is due to the 
fact that in these runs at lower resolution, the sum over 
P from the smallest wavenumber kmin = 1 to k\/2 in Eq. 
([5]) defines bands that are too narrow in Fourier space. In 
other words, it is linked to the lack of a well-defined iner- 
tial range in the simulations at lower Reynolds numbers, 
and we can only expect scaling to obtain in the limit of 
large Re- 

Figure 3] indicates that as the Reynolds number is in- 
creased, the contribution of the non-local interactions 
with the large scale flow to the total flux decreases 
(as well as the contribution of all non-local interac- 
tions, albeit at a slower rate). On dimensional grounds 
nLs(fcA) ^ Ulu\/L. Here, Ul is a characteristic velocity 
at the large scale L, and u\ is a characteristic velocity at 
the Taylor scale (note that this relation does not take into 
account that structures are in fact multiscale [l3|)- On 
the other hand, for Hls/H <C 1, we have Ii{k\) ~ u\/\. 

— 1/2 

As a result, we may expect I\.i,s{k\) /I{{k\) ^ Re . 
The condition Hls/H ^ 1 is not satisfied in the sim- 
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ulations at lower resolution, and it is unelear whether 
the departure of Run IV in Fig. [1] represents a conver- 
gence to a different scaling than ~ R^^'^ at very large 
Reynolds numbers. This will require further studies at 
higher numerical resolutions, a feat reachable with petas- 
cale computing. 

Other scaling laws can be observed in this series of 
runs; in particular, it is worth comparing the scaling of 
quantities for which data exist from laboratory experi- 
ments or from previous simulations. Figure [5] shows the 
Kolmogorov constant Ck as defined by the inertial range 
spectrum E{k) = Ck^^^^~^^^ ■ As a reference, we com- 
puted a best fit of the form Ck = clR\ + c, as suggested 
e.g. in [iOl, and obtained a = 4.60, h = —0.16, and 
c = 0.64. The value of c (that represents the asymp- 
totic value of the Kolmogorov constant for infinite Re) 
obtained from this fit is in good agreement with experi- 
mental results and atmospheric observations [iO, 13, al- 
though the values of a and h differ. We also note that 
the measured value of the Kolmogorov constant for the 
2048'^ runs is more than double the value of the expected 
asymptotic limit c, indicating that we are still far away 
from an asymptotic behavior for large R^. 

Figure [6] shows the skcwness 

S = {5uL{rf) / {5uL{TYf\ (9) 

and kurtosis 

K = {5uL{rf) / {5uL{rf)\ (10) 

of the longitudinal velocity increment = u • r/r 

(5Mi(r) = itL(x + r) - ul(x), (11) 

i.e. the component of the velocity in the direction of 
the increment. The skcwness and kurtosis were evalu- 
ated at two scales, r = A, the Taylor scale, and r ~ rj, 
the dissipation scale. In the latter case, only the results 
from runs III and IV show a dependence with R\ which 
is consistent with experimental results [i^. The behav- 
ior of these two runs further confirms that high Reynolds 
numbers are needed to observe scaling of turbulent quan- 
tities. 



IV. INTERMITTENCY AND STRUCTURES 

The Taylor-Green flows computed here correspond to 
an experimental configuration of two counter-rotating 
cylinders, studied in the laboratory for fluid turbulence 
as well as in the context of the generation of magnetic 
fields in liquid metals. These flows present both inhomo- 
geneities and anisotropics in the large scales, a resolved 
inertial range followed by a bottleneck, and a dissipative 
range. One may study the rate at which the symmetries 
of the Navier-Stokes equations are recovered in the small 
scales, and whether the statistical properties of the small 
scales are universal. In this section wc address the specific 
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FIG. 7: Scaling exponents using the ESS hypothesis in the 
1024=* and 2048^ runs. The scaling predicted by Kolmogorov 
and by the She-Leveque model are also given as a reference. 

question of the properties of the small scales through the 
evaluation of the anomalous exponents Cp of the longitu- 
dinal structure functions Sp of the velocity field, defined 
as: 

Sp ^ {6uL{r)P) - r^- , (12) 

assuming homogeneity and isotropy. In order to obtain 
better scaling laws, we use the Extended Self-Similarity 
hypothesis (ESS) [4^, |4^ in the particular context of 
plotting Sp as a function of S3. 

Figure [7] shows the scaling exponents (p in the 1024"^ 
and 2048'^ runs, computed using the ESS hypothesis. 
Similar results arc obtained without ESS and doing the fit 
only in the inertial range, defined as the range of scales 
where the so-called 4/5th law of Kolmogorov is satis- 
fied, namely 5*3 (r) ~ r. If we define stronger intermit- 
tency as stronger departure from the Kolmogorov scaling 
Cp = we note that as we increase the Reynolds num- 
ber, the intcrmittency increases as well, albeit slowly. 
Furthermore, for higher R\ (run IV), the departure from 
the She-Leveque model (47j increases (compared with run 
III), even for fixed values of p. The differences between 

for runs III and IV, albeit small, are at least one or- 
der of magnitude larger than the errors in the fit using 
ESS. As an example, in run III Ce = 1-746 ± 0.003 and 
Cs = 2.136 ± 0.007, while in run IV Ce = 1.7284 ± 0.0004 
and Cs = 2.0968 ± 0.0007. 

Here it is worth separating the discussion in two parts. 
On the one hand, the increase of the departure from the 
She-Leveque model as the Reynolds number and spatial 
resolution are increased indicates that the departure is 
not the result of lack of statistics. This change in the 
exponents for simulations with the same forcing at dif- 
ferent Reynolds numbers shows that huge Reynolds are 
required to obtain convergence of high order statistics. In 
fact, the larger the moment p examined, the larger the 
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FIG. 8: (Color online) Left: rendering of vorticity inten- 
sity in a small region of run IV. Only regions with > 
max{|a;|}/6.5 are shown (a; = V x v). Note the clustering 
of filaments into larger vorticity structures. The bars on the 
bottom indicate respectively the integral, Taylor, and dissipa- 
tion scales. Right: rendering of relative helicity in the same 
region (red is —1 and blue is 1). Only regions with absolute 
value larger than 0.92 are shown. 



relative difference between the C,p exponents measured in 
the two runs. On the other hand, it was shown in Rcf. 
that differences in the scaling exponents were mea- 
surable when considering two different forcings at simi- 
lar Reynolds numbers. These differences could be due to 
anisotropics in the flow, and in that case an SO (3) de- 
composition could be used to study whether the scaling 
exponents of the isotropic component of the flow are uni- 
versal. However, if there is a significant return to isotropy 
in the small scales, we then also expect the isotropic com- 
ponent to dominate when the Reynolds number is large 
enough. 

The intermittency of the flow is linked to the presence 
of strong spatially separated structures in the form of 
vortex filaments. The high R\ computation (run IV) dis- 
plays the same large-scale structure of bands as the run 
presented in 0. Conditional statistics analysis as the 
ones performed in Q keep showing a correlation between 
large scale shear and small scale gradients and enhanced 
intermittency. It has been noted by several authors that 
filaments tend to cluster into larger filamentary struc- 
tures; this is observed e.g. for supersonic turbulence 
and in the interstellar medium, and it has been ana- 
lyzed quantitatively in pioi] . When individual structures 
are studied in real space, filament-like clusters formed 
by smaller vortex filaments are observed here again (see 
Fig. [5]) , something that was not seen in simulations of the 
TG flow at lower resolution. This could be interpreted 
as a manifestation of self-similarity, and a more quanti- 
tative analysis will be presented elsewhere. In particu- 
lar, it would be of interest to compute the inter-cluster 
distance, and the intra-cluster intcr-fllamcnt distance, to 



see whether the space-filling factor of such flows diminish 
with increasing Reynolds number. Note that the vortex 
cluster reaches a global length comparable to the integral 
scale of the flow (indicated in Fig. [5]) ; as such, they may 
be a real-space manifestation of the trace of non-local in- 
teractions between small-scales (dominated by vortices) 
and large scales (dominated by the forcing) , giving a co- 
herence length to the flow. 

Figure [5] also shows the density of relative helicity 
v-a;(|v||a;|)~^ (w^ Vxv). Regions in blue and red cor- 
respond respectively to regions of maximum alignment or 
anti-alignment between the two fields (only regions with 
absolute relative helicity larger than 0.92 are shown). 
Note that regions with large relative helicity correspond 
to small vortex tubes, but the filament-like clusters have 
no coherent helicity. Regions with strong alignment fill 
a substantial portion of the subvolumc, even though the 
global (relative) helicity of the flow is close to zero. 



V. DISCUSSION AND CONCLUSION 

The data presented in this paper has allowed for a 
refined analysis of the behavior and structure of turbu- 
lent flows as the Reynolds number is increased. We have 
in particular showed that: (i) the bottleneck appears to 
have a constant width for the two higher i?e runs; hence, 
it is probably linked to the dissipation range, and to the 
depletion of nonlinearities as we approach this range; (ii) 
the scaling with i?e of the non-local energy fluxes, which 
indicates a weakening of non-local interactions as i?e in- 
creases. These first two results taken together point out 
to the fact that the bottleneck may not disappear in the 
limit of very high Reynolds number, since it has been ar- 
gued that its existence is linked to the relative scarcity of 
non-local interactions in Navier-Stokes turbulence, by op- 
position to, e.g., the magnetohydrodynamic (MHD) case. 
Indeed, when coupling the velocity to a magnetic field in 
the MHD limit, it was shown that the transfer of energy 
itself was non-local, and that the bottleneck was absent 
in numerical simulations of such flows; this can be under- 
stood in the following manner: as one approaches the dis- 
sipation range, few triadic interactions are available but 
in a flow for which the nonlinear transfer is nonlocal, the 
energy near the dissipative range can still be transfered 
efficiently to smaller scales since small-scale fluctuations 
arc transfered by the large scales [s^]. Finally, the de- 
parture of the anomalous exponents of velocity structure 
functions from standard models of intermittency such as 
the She-Lcvcque model seems to increase as the Reynolds 
number is increased. 

As noted before in [i^l , convergence to the asymptotic 
turbulence regime appears to be very slow: even though 
the nonlocal interactions do diminish with Reynolds 
number, they are still measurable at these resolutions. In 
run IV on a 2048^ grid at Rx ~ 1300, of the order of 10% 
of the energy flux is due to non-local interactions with the 
large scale flow, and the dependence of the energy flux 
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ratio Hls/II with Re for very large i?e is still unclear. 
This not only raises the question of the determination 
of higher order quantities at moderate Reynolds num- 
bers in simulations and experiments, but it also opens 
the door for a non-universal behavior of turbulent flows 
which may have to be studied in more detail than was 
previously hoped for. 
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